The Low-Field Critical End Point of the First Order Transition Line in YBa 2 Cu 3 07_5 
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We report on simulations of the first order phase transition in Y Ba-zCuzOt-t, using the Lawrence- 
Doniach model. We find that the magnetization discontinuity vanishes and the first order transition 
line ends at a critical end point for low magnetic fields in agreement with experiment. The transition 
is not associated with vortex lattice melting, but separates two vortex liquid states characterized by 
different degrees of short-range crystalline order and different length scales of correlations between 
vortices in different layers. 

PACS numbers: 74.20.De, 74.25.Dw, 74. 25. Ha 



Some features of the B — T phase diagram of clean un- 
twinned crystals of the high temperature superconductor 
YBa 2 Cu307_,5 (YBCO) are widely agreed upon. As the 
temperature is reduced the vortex liquid undergoes a first 
order phase transition to what is commonly assumed to 
be the Abrikosov vortex lattice state. There is strik- 
ing experimental evidence for this first order transition. 
Sharp drops in resistivity as well as discontinuities 
in magnetization and entropy occur simultaneously 
and mark the transition line. A neglected feature of the 
experimental first order transition line in YBCO is its 
termination at low (w 0.5T) magnetic fields. This has 
been consistently observed in all relevant experiments 
(i.e. the latent heat, magnetization jump and sharp re- 
sistance drop all disappear for fields smaller than some 
lower critical field) . Nevertheless the transition line 
is often shown extrapolated to T c , because a transition 
which separates phases of different symmetries, like a vor- 
tex lattice and a vortex liquid, cannot simply disappear. 

We report in this paper numerical results obtained for 
the Lawrence-Doniach (LD) || model for a clean super- 
conductor which show that the low field end point of 
the transition is not an artifact due to disorder (as com- 
monly thought) and that the interpretation of the first 
order transition as due to flux lattice melting must be 
incorrect. We see a first order liquid-liquid line with 
a critical end point at low fields, which in the B — T 
phase diagram shown in Fig.[j](a) agrees very well with 
the experimental YBCO "melting" line. On crossing the 
transition line, length and time scales increase discon- 
tinuously, but remain finite even in the low-temperature 
phase. The magnetization jumps we observe are shown 
in Fig.^(b) and are also in very good agreement with 
measurements in YBCO. Figure ^ is obtained using stan- 
dard YBCO values for the fitting parameters; viz for the 
Landau-Ginzburg parameter k=60, the mass anisotropy 
7=-\/m c /TO a h=7.5, the slope of the mean field transition 
line dB c2 /dT\ T =T c = -2T/K, the mean field T C =92.5K 
and the layer separation e?=11.4A. ft is perhaps note- 
worthy that one of us has frequently expressed doubts 
regarding the vortex-lattice melting scenario MM on the 



grounds that a finite temperature freezing to a flux lat- 
tice state is difficult to reconcile with the analytical result 
that off-diagonal long range order - a central feature of 
superconductivity - cannot coexist with a flux lattice at 
non-zero temperatures in three dimensions 0. 

The LD model for a layered superconductor consists of 
a stack of planes with Josephson coupling between neigh- 
boring layers. With the superconducting order parame- 
ter in the n th layer denoted as ip n , the Hamiltonian for 
the layered system with a magnetic field perpendicular 
to the layers is 
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FIG. 1. Phase diagram (a) and magnetization discontinu- 
ity (b) . Experimental data is taken from Ref. |2|Q (a) and 0] 
(b). N a t = vortices per layer (for numbers of layers see text). 
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where B = V x A, which we shall take as constant and 
uniform. In first approximation a(T) = a'(T — T c ) and 
02d(T) is constant; a',/?2D, J > 0. 

The model was simulated using Langevin (Model A) 
dynamics. The equation dtp/dt = —TdH/dip* + 77 was 
integrated numerically, using pseudo-random numbers to 
imitate the Gaussian white noise 77, (77* 77) = 2r k B T S r S t . 
We studied N a b vortices per layer for N c spherical lay- 
ers (see Ref. || for details of spherical boundary con- 
ditions), each of thickness do and spacing d and im- 
posed periodic boundary conditions ip\ = ipN c +i- For 
each layer ij) is expanded in eigenstates of the squared 
momentum operator (—ihV — 2eA) 2 . We retain eigen- 
states belonging only to the lowest eigenvalue 2eBh, 
(the lowest Landau level (LLL) approximation) which 
is a useful procedure over a large portion of the vor- 
tex liquid regime The model then depends on only 
two dimensionless parameters ]Tc| . The first one is the 
2D effective temperature and field parameter for each 
layer, given by arr = (doh/2e f3%r)B ksTY^aH, with 
«h = a(T) + eBh/rriab. The second, rj — J/\au\, con- 
trols the strength of the coupling between layers. 

If the order parameter varies only slowly across the 
layers |^| , the layered material behaves like a continuum 
with a mass anisotropy 7 and — P2D x d/do. The ef- 
fective mass m c in the c-direction can be expressed in 
terms of the layer coupling term: 77 = fi 2 / 2m c d 2 \aH \- 
In the continuum model of an anisotropic 3D type II su- 
perconductor, all thermodynamic properties of the sam- 
ple depend on only one parameter olt oc a B - With 
K = y^/3/2/j.Q x m a i,/eh we can express or in terms of 
measurable quantities as 



|a 2T r|| : 




FIG. 2. First order transition points in the «T-|«2Tf7| 
plane. N c varies between 8 and 80 for |q2t??| between 0.02 
and 2. The inset shows an example of the order parameter 
density p at constant |q2t7?|, here |a2T7?|=0.14, depending on 
increasing and decreasing olt- Squares and triangles corre- 
spond to 7V a i,=72 and jV a 6=144 respectively, filled symbols to 
cooling and open symbols to heating of the system. The tran- 
sition is marked by a clear hysteresis loop. Solid lines mark 
the discontinuity in p. 



_ ( V2h 3/2 Tr \ 2/3 [J_\ 2/3 dB c /dT\ Tc (T-T c )+B 
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Note that cct=0 corresponds to the mean-field H C 2(T) 
line and olt = — 00 to T=0. The LD Hamiltonian reduces 
to the continuum Landau-Ginzburg model in the limit 
77 — > 00, 0L2T — * 0, with a T — 77(2q!2t) 4 fixed. Because 
YBCO is closer to being 3D rather than 2D in character, 
we shall quote results in terms of olt rather than olit- 
The natural length scales of the continuum model are 
the magnetic length scale l m = \Jfij2eB^ which is pro- 
portional to the vortex separation distance in the layers, 
and the mean field coherence length £| | =%/ (2m c \ an \) 1 ^ 2 
parallel to the field. 

Instead of 77 we chose for our second independent pa- 
rameter the product la^T 7 ?! ■ Other than a factor 1 / VBT, 
\cxqtV\ contains only material constants and therefore 
varies slowly for rather a wide range of an and over con- 
siderable regions of the B-T plane. For comparison with 
experiment, we can express the coupling strength as 

I 1 ( *\ \ 1/2 1 1_ 

{a2Tm \8e 3 k B n ) K7 2 d3/2 ( BT )i/2- 

Note that the (unknown) 2D parameters of the layers in 
the model, K2D and do, cancel from the the definition of 
|cv2Tf?| if the relation k = y/d/do x K2D is used. Thus we 
only need to know the value k of bulk YBCO. 

Fig. |2| shows the phase diagram in terms of simula- 
tion parameters. As |a2T?7| increases and the system ap- 
proaches the continuum limit, the transition disappears 
at a critical point. This means that a first order tran- 
sition is not expected to occur in superconductors such 
as niobium for which the continuum approximation is 
appropriate. Our simulation results imply that the end 
point of the first order transition in YBCO is intrinsically 
not an effect of disorder. (Recent experiments [0 show 
that its location is significantly shifted to higher fields, 
(which corresponds to lower coupling |a2T^|) m twinned 
samples. This effect could be interpreted as an effective 
enhancement of the correlations along the field direction 
by correlated disorder. We expect point disorder to have 
only a slight effect on the location of the end point). 
Note that along the transition line ax is approximately 
constant, which means that the field and temperature 
dependence of the transition line behave like in a con- 
tinuum model where ax is the only scaling parameter in 
the system. 

The inset of Fig.|| shows an example of the kind of 
measurement used to locate the first order transition. 
The system displays hysteresis at it upon heating and 
cooling and the effects of this on the order parameter 
density are shown. The order parameter density is given 
by p — olt /3/2ira.H x (|V>| 2 ). The magnetization in the 
LLL model is 4irM = (/ioe 7 7./m a f ) )(|7/;| 2 ) (angular brack- 
ets signify a thermal average), which is in terms of our 
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simulation parameters 4%M — ir(B — B C 2(T))p/ oltk 2 , 
where B is the applied magnetic field. Thus we can work 
out the magnetization discontinuity from the discontinu- 
ity in p at the transition. The data points in Fig.|l](b) 
represent \a 2T v\ =1> 1-5, 2 and 2.5 at iV c =50, 60, 60/80, 
80. For these 4 transition points we have an average 
value of ax = —7.72, which yields the transition line in 
Fig.|l|(a). For |a 2 T?7|=2 and N c d/£\\ = 65 we see a clear 
transition, while for |a2T??|=2.5 and N c d/£\\ = 80 there 
is no sign of a transition in the range —8.3 < olt < —7.3. 

Near a critical end point we do not only expect the 
jump in the magnetization (order parameter density) to 
disappear, but we also expect there to be a divergence 
of the length scale of fluctuations in the order param- 
eter density of the system. We therefore looked at the 
density-density correlations of the order parameter near 
the critical point. A normalized density-density correla- 
tor is C d (Ar) = (|V>(r)| 2 |^(r + Ar)| 2 )/(|Vf) 2 - 1. Let 
us consider the case where Ar is a vector parallel to the 
c-axis. Plots of these correlations can be seen in Fig. [|. 
There is evidence of two length scales in the vicinity of 
the end-point. The short distance decay of the correla- 
tion function is dominated by the positional correlations 
of the vortices in the different layers. This length scale is 
mostly determined by olt and changes slowly in the vicin- 
ity of the critical point. However a second longer length 
scale becomes visible between olt = —7.6 and olt = —7.8 
as |oi2T??| is increased to its critical end point value. One 
can see in Fig. || that when |c«2T?/| =2.5 and olt — —7.8 
there is evidence of this much longer second length scale 
governing the decay of the correlation function at large 
distances. This length scale is associated with the density 
fluctuations at the critical end point and only become vis- 
ible once it is larger than that of the vortex correlations. 
Due to the small amplitude of these density fluctuations 
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FIG. 3. Density-density correlations along the c-axis near 
the end the first order transition line and the critical point 
(see Inset). Note the decreasing difference in correlations be- 
tween qt = —7.6 and «t = —7.8 as the transition disappears. 
For |q2t'7| — 2.5 and qt = —7.8 we see evidence for a long 
length scale associated with fluctuations in the average order 
parameter density. System sizes are N a i,— 72, 220 < N c < 260 
for a T = -7.6 and N c = 270 for a T = -7.8. 



very long simulation times are needed to determine the 
correlations within the statistical noise. 

It is often supposed that the first order transition in 
YBCO changes to second order below the end point, 
where no latent heat is visible but a "step" in the heat 
capacity C remains [| 1 1 1 . We believe however that this 
"step" can be identified with the onset of a small peak 
in the superconducting specific heat C s = C — C n , (n 
for normal state), which is known to arise from thermal 
fluctuations. This peak has been observed for example 
in niobium by Farrant and Gough jl^] . We find that the 
location and the height of the peak as well as the length 
of the rise (or width of the "step" ) in C from the low 
temperature value C S}Tn f (mf for mean field) to its max- 
imum agree well for the niobium and YBCO measure- 
ments taken from Ref . |Q a nd Q . The peak in C s in nio- 
bium obeys LLL scaling jl2j] and is found at olt /Pa = — 6, 
where 0a = 1.16, i.e. ax ~ —7. The data which shows 
the specific heat in YBCO is given as C minus C(B = 0). 
The latter is near the "step" approximately equal to the 
low temperature value C s ^ m f + C n , so that the plotted 
quantity is approximately C s — C SjI „/. The maximum oc- 
curs for example for B = 0.25T at T « 91. 4K which cor- 
responds to ax — —7.2. The width of the "step" in nio- 
bium Aa T ~ 2. In YBCO for B = 0.25T the specific heat 
rise associated with the step takes place in the tempera- 
ture region 91 — 91.4A", which corresponds to Aolt = 2.8. 
For niobium C s is at its maximum 5% larger than C s . m f. 
In YBCO we have to divide the plotted data by C s , m / to 
compare with this value. C S;m / is roughly given by the 
step in C at the zero field transition which we take from 
Ref. g. We find that (C s - C s , mS )/C s , m} » 0.02, which 
is of the same order as in niobium. The specific heat 
"step" in YBCO at different fields has approximately the 
same amplitude as well as width and position when ex- 
pressed in terms of olt-, i.e. the "step" feature obeys LLL 
scaling. The quantitative agreement between YBCO and 
niobium strongly suggests that we are dealing with the 
same phenomenon and therefore that there really is no 
sharp specific heat step in YBCO. 

The existence of the critical point implies that no sym- 
metries are broken at the transition, which means it can- 
not be a liquid-crystal transition. We find indeed that the 
vortex matter is liquid on both sides of the transition. In 
Fig. ^ we show examples of density-density correlations 
above and below the transition for |cv2T?7| = 1 which cor- 
responds to a transition temperature of 83 K in YBCO. 
Fig. ^(a) shows the correlations along the c-axis as previ- 
ously seen in Fig. |^. We see an exponential decay of the 
correlation function with a finite length scale l c for den- 
sity correlations below as well as above the transition. 
Only for a liquid phase would this correlation function 
have an exponential decay. In Fig. ^(b) we show exam- 
ples of measurements of the 2D Fourier transform of the 
density-density correlator Cd for Ar parallel to the lay- 
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ers, normalized by its high temperature limit, which is 
essentially the structure factor of the system (see Ref. || 
for details). It has its first and largest peak at the first 
reciprocal lattice vector of a triangular lattice, k w 2.6 in 
units of magnetic length. The inverse width at half max- 
imum of the Lorentzian fits to this peak gives the length 
scale of crystalline order in the layers, l a b- The length 
scales taken from fits in figure ^ are in order of decreas- 
ing a T l c /i\\ =13, 15, 20, 28, 46 and l ab /l m =1.33, 1.57, 
1.94, 2.47, 2.73. 
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FIG. 4. (a) Density-density correlations along the c-axis 
(with linear fits), (b) structure factor (with Lorentzian fits) 
at ja 2 T7?j = l. The system size is N ab =72, 7V C =270. 

The length scales have a discontinuity at the transi- 
tion. This jump is found to grow with distance from 
the end point, as one would naturally expect. We also 
see a rapid growth of the length scales below the tran- 
sition. An exponential growth of length scales with 
| ar | 3 / 2 ~ (T c — T) 3 / 2 , has been predicted for the low tem- 
perature regime from perturbative studies around zero 
temperature ||. We can extrapolate the growth in the 
crystalline length scale obtained from Fig. |] (b) just 
below the transition assuming exponential growth with 
|ax| 3 / 2 - For a decrease of Aax ~ 1-2, which corresponds 
in YBCO to cooling by only IK below the transition at 
83K, we obtain an increase in the range of crystalline or- 
der by a factor of 3. It is therefore likely that not far 
below the first order liquid-liquid transition length scales 
reach the system size in a pure system or in real crystals 
a "Larkin" like length scale (dependent on the amount 
of disorder present) so that the vortex liquid is corre- 
lated and effectively crystalline over large length scales. 
Although the structure factor in a liquid is rotationally 
symmetric, coupling with the underlying lattice may for 
long length scales lead to the appearance of Bragg-like 
peaks 1 1 3 1 - 

The longest time scales in the system, given by the 
decay of 3D Fourier component of Cd at the first recip- 
rocal lattice vector in the a6-plane and fc=0 along the 
c-axis (not shown), increase discontinuously across the 
transition. Such behavior may explain the sharp fea- 
tures of transport coefficients like resistivity. We also see 
extremely fast further growth of time scales below the 



transition in agreement with the fast decay of resistivity 
to zero as the temperature is lowered JTo| ] . 

From previous simulations using the same model with 
periodic boundary conditions (PBC) in all directions in- 
stead of the geometry of spherical layers in a radial field, 
a vortex lattice melting transition is reported Jl4],[l5| . We 
find disagreement in the location of the transition due to 
the different choice of boundary conditions only for low 
|o!2T?7|- We believe that our choice of boundary condi- 
tions is more likely to reflect the real physics |l(| . In the 
3D-like regime appropriate to YBCO, the LLL-LD model 
may well show the same behavior for PBC as it does for 
spherical layers. For couplings high enough to see the 
critical point, the LLL-LD model has to our knowledge 
never been investigated using PBC. The largest system 
sizes used in the simulations with PBC are of the order 
of 40 vortices x 20 layers, smaller than the ranges of cor- 
relations we find below the transition. This would make 
the vortex liquid indistinguishable from a vortex lattice. 

We have also attempted to compare our simulations 
with experimental data on Bi2Sr2CaCu20s, for which 
the first order transition is in a region of the phase di- 
agram where the LLL approximation is not good. As 
might be expected our results are not as quantitative as 
they are for YBCO but they are none the less still useful. 
Details will be given in [To[ . 

In summary, we have found in a simulation the low 
field end point of the first order transition line in YBCO 
and obtained results in excellent agreement with exper- 
iment. The existence of a critical point implies that the 
vortex matter is liquid above and below the transition, 
and we were able to observe this directly in our simula- 
tion. Our results suggest that the transition in YBCO, 
which is commonly interpreted as vortex lattice melting, 
is of a liquid-liquid nature and that the vortex crystal 
state does not exist. 
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